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Q,' ABSTRACT 

< 

^ ■ The 7-ray burst (GRB) 060218/SN 2006aj is a peculiar event, with the sec- 

ond lowest redshift, low luminosity, long duration, chromatic lightcurve features, 
^ ■ and in particular, the presence of a thermal component in the X-ray and UV- 

'^ ■ optical spectra. Thanks to detailed temporal and spectral coverage of the Swift 

lO ' observatory, the abundant data allow the GRB prompt emission to be modelled 

in great detail for the first time. The low flux of prompt UV/optical emission 



^ . disfavors the conventional internal shock/synchrotron radiation models, which 

->^ . generally predict strong UV/optical emission. Here we show that the unusual 

p . prompt emission of GRB 060218 can be produced by inverse-Compton scattering 

Q . of shock-accelerated relativistic electrons off the detected thermal photons. A 

^ . pair of (forward plus reverse) shocks form when a relativistic outflow interacts 

c^ . with a preexisting slower shell. The observed 7-ray emission and X-ray emission 
arise from the reverse-shocked and forward- shocked regions, respectively. A fit 

K^ , to the data requires an initially increasing outfiow luminosity, which is consistent 

^ ! with the prediction of the popular collapsar model. 

Subject headings: gamma-rays: bursts — relativity — shock waves 



1. Introduction 

The 7-ray burst (GRB) 060218 was detected by the Burst Alert Telescope (BAT) on- 
board Swft satellite on 2006 February 18.149 UT (Cusumano et al. 2006). The prompt 7-ray 
emission lasted about 2000 s and thus it is one of the longest GRBs. Spectroscopic measure- 
ments after the burst revealed that this event was associated with a type-Ic supernova, SN 
2006aj (Modjaz et al. 2006; Plan et al. 2006; Sollerman et al. 2006; Mirabal et al. 2006; 
Cobb et al. 2006) and that its redshift is z = 0.0331 (Mirabal & Halpern 2006). This is the 



-2- 



second lowest redshift GRB after GRB 980425/SN 1998bw. The burst has several unusual 
7-ray and X-ray properties (Campana et al. 2006, see Fig.l). First, both prompt 7-ray emis- 
sion and X-ray emission detected by BAT and the X-ray Telescope (XRT) exhibit a slow 
brightening and a final steep decay, unlike spiky light curves of classical bursts at higher red- 
shifts. The lightcurves are chromatic, with the gamma-rays fading up earlier than X-rays. 
Second, the peak 7-ray fiux is as low as ~ 10~^ erg cm~^ s^^ and the isotropic-equivalent 
7-ray energy is a few times 10^^ ergs, about three orders of magnitude lower than the typical 
energy of classical GRBs. Third, a soft thermal component was detected by XRT ~ 300 
to 10^ s after the burst, which a luminosity ~ 10^^ erg s^^. The thermal component is also 
identified in the UVOT data at ~ 10^ s. This soft component is likely due to the break- 
out of a radiation-dominated shock, possibly driven by a cocoon (Meszaros & Rees 2001; 
Ramirez- Ruiz et al. 2002; Zhang et al. 2003), parts of uncorked envelope matter (Waxman 
& Meszaros 2003) or parts of accelerated envelope matter (Colgate 1974; Tan et al. 2001). 
Fourth, the prompt optical/UV emission is extremely faint and its fiux is about five orders 
of magnitude fainter than that of the prompt X-ray emission. Finally, it is interesting to 
note that the observed X-ray afterglow decayed as a simple power law with time (~ t^^'^) 
up to 10 days after the burst. 

In the conventional internal shock models (for reviews see Zhang & Meszaros 2004 and 
Piran 2005), collisions among relativistic shells with different Lorentz factors would lead 
to internal shocks and the prompt emission of a GRB could arise from shock-accelerated 
electrons by synchrotron radiation and/or inverse-Compton scattering (e.g. Rees & Meszaros 
1994; Dai & Lu 2002; Zhang & Meszaros 2002). In these models, the typical synchrotron 
self-absorption frequency is in the optical/UV band, so the prompt optical/UV emission 
fiux is about three orders of magnitude fainter than the prompt X-ray fiux because the 
synchrotron spectrum is uF^ oc u^^^ in the optical/UV to X-ray bands. The models predict 
a strong prompt optical/UV fiash (Meszaros & Rees 1999), especially in low-redshift bursts. 
The tight constraints of the prompt UV-optical data of GRB 060218 therefore implies that 
the conventional internal shocks/synchrotron radiation models are inapplicable at least for 
this burst. 

In this Letter, we propose a model in which a relativistic outfiow interacts with a 
preexisting slower shell, leading to a pair of relativistic (forward plus reverse) shocks. The 
slower shell could have been a part of the stellar envelope ejected from the central engine 
prior to the GRB or a circum-burst dense clumpy cloud. We show that the unusual prompt 
emission of GRB 060218 could be produced by inverse-Compton (IC) scattering of shock- 
accelerated electrons off the detected shock breakout thermal photons. To our knowledge, 
this is the first detailed theoretical model of the prompt emission of a GRB, which is made 
possible thanks to the broad temporal and spectral coverage of GRB 060218 by the Swift 
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satellite. In §2 we present the IC spectrum of a relativistic shock in a general case. In §3 
we discuss the shock dynamics and calculate the emission luminosity and in §4 we constrain 
the model parameters for GRB 060218. Our results are summarized in §5. 



2. IC spectrum in relativistic shocks 

We first discuss a general case in which a relativistic shock with bulk Lorentz factor of 
7 occurs simultaneously with a soft photon source. We assume that the luminosity of this 
photon source is Lq and its radius is much smaller than the shock radius R. We now consider 
IC scattering of the shock-accelerated electrons with the soft photons. If the shock expands 
outwards isotropically, there should be many head-on scatterings between the electrons and 
the soft photons so that some electrons would quickly cool down, as argued by Brunetti 
(2000), Wang, Li & Meszaros (2006), Wang & Meszaros (2006), and Fan & Piran (2006). 
Letting eg be the typical energy of the soft photons in the observer's frame, we have the IC 
power and characteristic energy from a relativistic electron with Lorentz factor 7e, 

Pile) = arc^VeK, (1) 

e,(7e)=7e%- (2) 

The factor 7^ in eq. (1) is introduced to transform the power in the shock rest frame to the 
observer's frame, ay is the Thomson cross section, and Uq = Lo/(47ri?^7^c) is the energy 
density of the soft photons in the shock rest frame. If eq. (1) is multiplied by a factor of 4/3, 
one obtains the IC power of a relativistic electron in an isotropic photon field. In the present 
case, a beam of soft photons illuminates the shock-accelerated electrons along the radial 
direction. One needs to introduce a correction factor to the IC power in an isotropic photon 
field (Fan & Piran 2006). This correction combined with the factor 4/3 gives a coefficient 
of order unity, which is neglected in eq. (1). The peak spectral density (i.e., power per unit 
energy) is approximated by 

_ Pile) _ (Tt Lq 

e^(7e) AttR^ eo 

which is independent on 7e and (apparently) on 7. This expression represents the number 
of scatterings with soft photons per unit time for a certain electron in the Thomson limit. 
Thus, the peak spectral luminosity is given by 

where Np is the total number of electrons. 
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Similar to Sari, Piran & Narayan (1998), if a relativistic electron loses its energy through 
the IC discussed here, the critical Lorentz factor 7c above which electrons cool in the ob- 
server's time t is defined through 77cmeC^ = P{lc)t, i.e. 

rUeC 47rmeC^7-R^ ,_, 

Tc = — 7777 = ri — • (^) 

As usual, we assume that the instantaneous shock-accelerated electron energy distribution 
follows a power law: dNejAle oc 7~^ for 7e > 7m. If a constant fraction (eg) of the post-shock 
energy density goes to electrons, one has 7^ = (^p/'"^e)ee[(p~2)/(p — 1)]7 = 610ee5'p7, where 
g, = 3{p-2)/{p-l). 

If the cooling time of electrons with Lorentz factor of 7^ is shorter than t (i.e. 7^ > 7c), 
the steady electron energy distribution turns out to be dNe/d'^e oc 7"^ for 7c < 7e < 7m and 
dNf./d'^e oc 'je for 7e > 7m. In such a fast cooling regime, the spectral luminosity at e^ 

is (e.g. Blumenthal & Gould 1970) 

T J -^e-Y,maxl,f7/fcJ j ^c ^ ^ 7 ^ Cm,; /^\ 

^^ ~ I L (e U Y^l'^ie le Y'p'^ e > e ^ ' 

(_ '^ej,nia.x\'^m/ '^cj V 7/ m,/ 5 ^^7 _ ^my 

where ec ~ 7^60 and em ^ 7meo. 

On the other hand, if 7^ < 7c (slow cooling), the steady electron energy distribution 
becomes dN^/d'^e oc 7~^ for 7m < 7e < 7c and dN^./d'^e oc 7,7 for 7e > 7c. In this case, 

the spectral luminosity at e^ reads 



-'^e-y,maxV"-7/ "-cj j '^m ^ C'y ^ Cc 

r ff /f ^-(p-i)/2('f /^ ^-p/2 f > f 



(7) 



3. Shock dynamics and emission luminosity 

We assume that the luminosity of a relativistic outfiow is Lw{t) oc t^ for t < t^, 
and L^{t) oc t~^/^ for t > t^, as is suggested by the popular collapsar model for long- 
duration, soft-spectrum GRBs (MacFadyen, Woosley & Heger 2001). The Lorentz factor 
of the outfiow is 7^ ^ 1. The co moving proton number density in the outfiow rest frame 
is n^ = Lw / {ATiR^'-yf^nipC') . For simplification, we also assume that a preexisting shell is 
so slow (as compared with the outfiow) that the shell can be considered as being at rest, 
and its uniform proton number density is n. We assume that the collision happens at an 
initial radius Ri. Similar to Sari & Piran (1995), the outfiow-shell interaction beyond Ri 
could be described through two shocks, a reverse shock that propagates into the outfiow. 



-5- 



and a forward shock that propagates into the shelL There are four regions separated by the 
two shocks, (1) the unshocked shell, (2) the forward- shocked shell, (3) the reverse-shocked 
outflow, and (4) the unshocked outflow. The shock dynamics depends on whether the re- 
verse shock has crossed the outflow. It is clear that at two stages, i.e., during and after 
the reverse-shock crossing, the shock dynamics is different and thus the emission luminosity 
should also be different. In the following, we discuss these two stages separately. 



3.1. During the reverse-shock crossing 

As shown by Sari & Piran (1995), if 7^ 3> ^«,/n, the reverse shock is relativistic and 
the Lorentz factor of regions 2 & 3 measured in the rest frame of region 1 is 



/7^\i/2 /n^\i/4 1/8 -1/8., ,, 



where the convention Q = 10^ x Q^ in cgs units has been adopted and tj = tj 3 x 10^ s = 
Ri/{2-flc). The shock radius becomes i? = 2-flct+Ri = 2-flc{t+ti) = 6.3xlO^^L^'^Jo<^^^(^3+ 
^j.s)^ cm. In addition, the Lorentz factor of region 3 measured in the rest frame of region 4 
is 

"'" = (t)'" i^y" = l-54i;'/„'"J"7»,i(*» + *«)■'*. (9) 

From eq. (5), we calculate the cooling Lorentz factor of the relativistic electrons 

7e = 2.0L^^,,Lf,,nf/%\t, + U,s?^\ (10) 

The characteristic Lorentz factors of the relativistic electrons in regions 2 & 3 are, respec- 
tively, 

20ee-29pLli^5ong^^^it5 + ti,3y^^^, region 2, 
^O^e-29pLZ^ionJ^lw,i{t3 + ti^sY^^, region 3. 

for a soft photon source with a black body spectrum and a typical photon energy cq ~ ksT ~ 
0.16keV (as in GRB 060218), the cooling energy of the scattered photons is 

ee = OML^',,Lf,,n-'\Ttf{ts + U,,f'^ keV, (12) 

where Cr = /csT/O.lGkeV. The characteristic energies of the scattered photons from regions 
2 & 3 are, respectively, 

Uel_^glL%n^'^'\T{h + t.s)"^/' keV, region 2, ^^^^ 

16ee,-2^p^^,5'o^^9^S^,iCT(t3 + h^f/^ keV, region 3. 



Im — \ ,„ :r-l/8 1/8 /, . , M/^ . „ \^^J 
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The total electron numbers in regions 2 & 3 are A^e,2 = 47ri?^(27|,ct)n = 3.1 x lO^'^LjgQng t3(t3+ 
tj,3)^/^ and A^e,3 = A^e,2(7rs/7fs) = 1-5 X 10^^L^'^5on9''%^,it3(t3 +ti,3), respectively. From eq. 
(4), we then obtain the peak spectral luminosities for regions 2 & 3, 

2.7 X lQ*^Lo^,^L%nl'\;^%{h + U^,)-^'^ ergkeV-^ s-\ region 2, ^^^^ 
1.3 X 10^^Lo,46n97^^iCr^^3 ergkeV^"*^ s^"*^, region 3. 

According to equations (12)-(14), the spectral luminosities at scattered photon energy e^ in 
regions 2 & 3 are 

^ 2.2xl046L:/|on9^V^'(^3 + t,3)^/VikevergkeV-is-i, region 2, 

^^ ^ 1.0xl0^^4Vf7.,iCT'^'(^3 + t,3)^/\-;/k'evergkeV-is-i, region 3, ^ ^ 

for ec < e^ < em, and 

^ r 4.9 X 10^^er-^2^r^i^S^^'n9"^'~'^^'C?''^/'(t3 + h^)-^'-'^'%1iy ergkeV-^ s-S region 2, 
^^ \ 7.5 X 10^6er_^2^r^4y/^?^'^/'7:,iC?-'^/'(t3 + t,3)(^+^)/^e;f V ergkeV^ s-\ region 3, 

(16) 
for ec < em < e-y. Here e^^ikev = e^/lkeV, and the coefficients in eq. (16) have been 
calculated for p = 2.5. 



3.2. After the reverse-shock crossing 

We assume that the reverse-shock crossing time is t^ ~ tj. After this time, the central 
engine luminosity is insignificant because it decays as L^ oc t~^^^ (MacFadyen et al. 2001). 
Although the electrons in region 3 have cooled down, the higher-latitude scattered photons 
would reach the observer at a later time. This is the well-known curvature effect (Kumar & 
Panaitescu 2000; Zhang et al. 2006; Wu et al. 2006; Liang et al. 2006). In the relativistic 
limit, this effect leads to the emission luminosity L^^ oc t^(2+/3) (^-^j^ere P is the spectral 
index). As argued by Zhang et al. (2006), the temporal index a = 2 + j3 conclusion is 
essentially unchanged by deceleration of region 3. Meanwhile, the forward shock still sweeps 
up its surrounding shell and its Lorentz factor decreases as 7fs ex: t~^^^. The cooling Lorentz 
factor decays as 7c oc (7fsf/Qt)~^ oc 7fs-R^t^^ oc t~'^/® (see eq. [5]) and the cooling energy 
of the scattered photons deccay as ec oc t^^/"^. The peak spectral luminosity increases as 
Le^,ma.x OC NeR'"^ OC i? OC t^^^. Thus the spcctral luminosity at ec < e-y < em decays as 

i^e =^.„max(e^Ac)~'/'0Cr5/8. 
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4. Application to GRB 060218 

The prompt emission liglit curves detected by BAT and XRT are presented in Campana 
et al. (2006): Tlie 15-150 keV luminosity and the 0.3-10 keV luminosity ~ 900 s after 
the burst increases as Lbat oc t^T-^ with a^^i ~ 0.3 and Lxrt c<: t"^-^ with a^^i ~ 0.7, 
respectively. Subsequently the 15-150 keV luminosity decays as Lrat c<: t"^'^ with 0:^2 ~ 
—4.5, while the 0.3-10 keV luminosity declines slowly as Lxrt oc t"""'^ with ax,2 with ax,2 ~ 
—0.6 until t ~ 2600 s, after which the lightcurve declines rapidly as Lxrt c< t"""-^ with 
(^x,3 ~ —4.0. The different light curves detected by BAT and XRT suggest that their 
corresponding emission could have originated from different regions. In the following, we 
show that the X-ray emission and 7-ray emission would be produced from regions 2 and 
3 in our model, respectively. An example fit to the observational data with our model is 
displayed in Figure 1. 

We first consider the emission luminosity at the peak time tpeak ~ 900 s. Taking p ~ 2.5, 
ti ~ 10^ s and ksT ~ 0.16 keV, we calculate the 15 — 150 keV luminosity by using the second 
expression of eq. (16), 



/.150keV 

Lbat = / L,-,de^ ~ 3.3 x lO'\l-',,L''Jf,nf'^l',eTgs-' ~ 1.8 x lO^^rgs'^ (17) 

Jl5keV 



"150 keV 
'15keV 

Using the first expression of eq. (15), we have the 0.3 — 10 keV luminosity 

/•lOkeV 

Lxrt = / L,-,de^ ~ 1.7 x 10^' L^J^^nl^'' ergs'' ~ 3.3 x lO^^ergs-^ (18) 

^0.3 keV 

From eqs. (17) and (18), we find the shell number density and the outflow luminosity at the 
peak time tpeak - ty,, 

ng ~ 0.02e-i275, (19) 

L»,50 ~ 0.3el/J_^-fi;i, (20) 

respectively. 

Second, we discuss the light curves before tpeak- Theoretically, the 7-ray and X-ray 
luminosities increase as Lrat oc Lw and Lxrt oc LJ , respectively. Considering L^ oc t'^ 

before t^, we require a^^i ^ (6 — p)k/8 ~ 0.3 and a^^i — 7k/8 ~ 0.7, giving k ~ 0.8. This 
indicates that the outflow luminosity increases slowly with time before tpeak- In the collapsar 
model, numerical simulations by MacFadyen et al. (2001) indicate that the accretion rate 
of the fallback disk by the central stellar-mass black hole typically shows a slowly-rising 
behavior before the well-known M oc t^^/^ accretion starts. The outflow luminosity has a 
similar temporal power because of a nearly constant Blandford-Znajek power efficiency for 



the rotating black hole (McKinney 2005). Thus, our fit is consistent with the prediction of 
the collapsar model. 

Third, after t^, the reverse shock has crossed the outflow and the curvature effect of 
region 3 leads to the sharp 7-ray luminosity decay, i.e. Lbat oc t^(^+'^) ~ t^^-^^. This is 
marginally consistent with the steep decay of the 7-ray emission. A larger value of J9 or a 
potential spectral break at higher energies would make the curvature effect more consistent 
with the observed decay. Meanwhile, region 2 could have entered the self-similar solution 
of Blandford & McKee (1976) and its X-ray luminosity decays as Lxrt oc t"^'^. This is in 
good agreement with the observed value of ax,2- 

Fourth, we note from Campana et al. (2006) that after ~ 2600 s both the non-thermal 
component luminosity and the thermal component luminosity decayed rapidly. This steep 
decay is the result of the combination of the rapid decay of the cocoon emission (Pe'er et al. 
2006) and the curvature effect of region 2 after the forward shock has crossed the shell. 

Finally, the forward shock, after crossing the slow shell, would sweep into the circum- 
burst stellar wind and produce a normal afterglow as observed by Swift XRT. A single power 
law light curve of this afterglow indicates that the outflow is not strongly collimated (Fan, 
Piran & Xu 2006). 



5. Conclusions 

The abundant multiwavelength prompt emission data of GRB 060218/SN 2006aj make it 
possible to constrain the prompt emission mechanisms and model parameters in great detail. 
The data clearly disfavor the conventional internal shock/synchrotron radiation models that 
predict strong prompt UV optical emission. Instead the data are consistent with a model 
that invokes inverse Compton scattering off the thermal photons due to shock breakout. 
The chromatic lightcurve features of the prompt emission require that the prompt BAT- 
band emission and XRT-band emission are produced in two different regions, and we identify 
these as a pair of shocks upon the interaction between a relativistic outflow and a preexisting 
slow shell. Fitting the data with model poses constraints on the physical parameters of the 
outflow and the shell, in particular, requires a central engine wind with a slowly-increasing 
luminosity. This behavior is consistent with the prediction of the popular collapsar model. 

This work was supported by NASA under grants NNG05GB67G, NNG05GH92G and 
NNG05GH91G (BZ, ZGD & EWL) and the National Natural Science Foundation of China 
under grants 10221001 & 10233010 (ZGD) and 10463001 (EWL). 
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Fig. 1. — An example fit to the prompt emission data of GRB 060218 witli our model, for 
p ~ 2.5, k ~ 0.8, and t^ ~ 900 s. The red dots and green squares represent the XRT and BAT 
data (taken from version 2 of Campana et al. 2006), respectively. In our fit, the required peak 
outflow luminosity is L^ 50 ~ 0.3eJ_2^J^, the shell number density is rig ~ 0.02e7!^27ioi5 ^^d 
the shell width is Ai?shcii ~ 3.2 x lO^^Cg _27^ ( cm. The rising segments of our theoretical 
7-ray and X-ray light curves are Lbat oc t^^^'P")^/^ and Lxrt oc t^^^^^ respectively. During the 
period between ~ 900 to 2600 s, the decaying index of the X-ray luminosity is —5/8 and 
the decaying index of the 7-ray luminosity becomes — (2 + /3) because of the curvature effect 
after the reverse shock has crossed the outflow. At ~ 2600 s, the forward shock has crossed 
the shell, and hence, the X-ray luminosity shows a rapid decay due to the curvature effect. 
However, the actual decay becomes steeper because of the rapid decrease of the observed 
thermal luminosity. 



